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Abstract 


The most useful characterization of a gravitational wave detector's 
performance is the accuracy with which astrophysical parameters of 
potential gravitational wave sources can be estimated. One of the most 
important source types for the Laser Interferometer Space Antenna (LISA) 
is inspiraling binaries of black holes. LISA can measure mass and spin to 
better than 1% for a wide range of masses, even out to high redshifts. The 
most difficult parameter to estimate accurately is almost always luminosity 
distance. Nonetheless, LISA can measure luminosity distance of 
intermediate-mass black hole binary systems (total mass~10 4 M@) out to 
z~10 with distance accuracies approaching 25% in many cases. With this 
performance, LISA will be able to follow the merger history of black holes 
from the earliest mergers of proto-galaxies to the present. LISA’s 
performance as a function of mass from 1 to 10 7 and of redshift out to 
z~30 will be described. The re-formulation of LISA’s science requirements 
based on an instrument sensitivity model and parameter estimation will be 
described. 



The Problem 


■ Gravitational wave detection is “capability-driven, ” meaning 
that the gravitational wave science that we want to do is 
determined by the measurement capability that we can 

muster. 

■ However, circumstances compel us to talk about the science 
that we can perform to justify. 

■ NASA requires that the instrument performance requirements 
be strictly derived from the science requirements. 

• How it actually happened: 

- Conceive an instrument that seems promising 

- Assess what the instrument can do 

- Figure out what sources can be observed with that capability 

- Figure out what science one can do with those observations 

- Iterate, endlessly. 

■ So, what can one do with a GW detector, such as USA ? 



The LISA Sources 


Ultra-compact binaries 
• ~1M & 

• Galactic and extragalactic 

• 1000’s - 10,000 

• Confusion foreground 



Extreme mass-ratio 
inspirals 
- 10 / 10 6 M & 

z < 1 

10’s - 100 per year 


Massive and 
intermediate-mass black 
hole binaries 

• 10 2 - 10 7 M @ 

• z <20 

• 10’s to 100 per year 


Cosmological back- 
grounds, bursts and 
unforeseen sources 
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GW Observations 


■ Concentrate on inspirals, rather than burst, quasi-periodic, or 
stochastic sources. 

■ With few exceptions, we don’t know the sources beforehand 

- We only know the general form (template) of the signals to look for. 

- The parameters of specific sources are unknown prior to observation , and other 
characteristics , such as the event rate, are poorly known. 

- The properties of the sources are determined in the process of “detection” 
through template matching. 

■ LISA can determine source parameters pretty well (by astrophysics 
standards), especially quantities encoded in frequency. 

■ To a significant degree, we’re going to discover what science can be 
done with our GW detection as we go along. 

■ So, how do we characterize the discovery space? 

• How do we show what sources can be observed, should they exist? 






“Detection” vs. “Observation” 


• “Detection” means identifying the presence of a source 

- Determine the likelihood of a source signal in the data 

- Figure of merit is SNR 

- Detection can be useful, but doesn’t answer specific scientific 
questions, like what is the merger rate of intermediate mass black 
holes at z- 10? 

■ “Observation” means getting useful estimates of source 
parameters 

- Estimate the value of source parameters 

- Figure of merit is uncertainty in parameter estimation, e.g., luminosity 
distance, masses, spin parameters, sky location, etc. 

- Observation is really the desired metric of scientific capability. 



Detection with LISA 
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Observation with LISA 


■ The range of sources is large. 

■ Our science objectives are diverse. 

■ How do we know that the instrument we’ve designed can 
satisfy our science objectives? 

■ What other science objectives might it satisfy? What’s the 
discovery space? 

■ Here’s the process: 

- Write down the objectives 

- Formulate the desired science observations (which we think that we 
can fulfill) which will satisfy the objectives, with requirements on the 
estimation of source parameters 

- Calculate a performance model of the instrument 

- Analyze whether that instrument can perform those observations. 

• And accommodate the possible range of source/detector 
interactions. 



The Unknown Source 


■ Waveform depends on 17 source parameters 

- Intrinsic parameters of the source: m p m 2 , s p s 2 , eccentricity, separation, 
orbital phase. [Chirp mass and reduced mass.] 

- Extrinsic parameters of the observer: luminosity distance, sky position, 
orientation of the source orbit on the sky, orientation of the detector, merger 
time, orbital phase of the detector. 
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Evolution of the characteristic strain signal 



log obtferved frequency [H*] 


Characteristic strain for various masses (1 0 5-7 , 10 3 M 0 ) and redshifts (1 ,3,5,7), 
1 yr, 1 mo and 1 day. (Sesana, Haardt, Madau and Volonteri, 2004 ) 





Instrument Assumptions 


■ Starting frequency for integration 

- Lower edge of detector sensitivity (3x1 0 5 or 1x1 0 4 Hz) 

- Time before coalescence 

■ Number of interferometers 

- 6 links (all 3 arms) 

- 4 links (1 equivalent Michelson) 









Science Objectives and Investigations 


■ LISA science spans seven science objectives 


Understand the formation of massive black 
holes 

Confront General Relativity with 
observations 

Trace the growth and merger history of 
massive black and their host galaxies 

Probe new physics and cosmology with 
gravitational waves 

Explore stellar populations and dynamics in 
galactic nuclei 

Search for unforeseen sources of 
gravitational waves 

Survey compact stellar-mass binaries and 
study the structure of the Galaxy 



■ Each science objective has one or more science investigations by 
which it will be achieved. For example, 


Trace the growth and merger history of 
massive black holes and their host galaxies 

Determine the relative importance of 
different black hole growth mechanisms as 
a function of redshi ft 


Determine the merger history of IxlO 4 to 
3x10 I s M & black holes from the dawn of 
galaxies (z~20) to the era of the earliest 
known quasars (z~6) 


Determine the merger history of 3x1 0 s to 
IxlO 7 M @ black holes at later epochs (z<6) 




Science Requirements 


Science Investigation 


4.2.1 Determine the relative importance of different 
black hole growth mechanisms as a function of 
redshift 


Observation Requirement 



OR2.1: LISA shall ha 1 
massive black hole bin 
mass in the range 3x1 C 
smaller mass in the rar 
z = 10, with fractional 
25% for luminosity dis 
for spin parameter at n 
maintain this detection 
increase the number oi 

ve the capability to detect 
lary mergers, with the larger 
K M @ <Mj< 3x1 0 5 M & and a 
lge 10 3 M @ < M 2 < 10 4 M & at 
parameter uncertainties of 
stance, 10% for mass and 10% 
raximal spin. LISA shall 
i capability for five years to 
f observed events. 



0 ; ] 00 1 00 
l i ■ !d>: vs) 



Frequency (Hz) 



Parameter Uncertainties 

— i 

M, 

m 2 

D l Uncertainty 

Spin Uncertainty 

SNR 

1.00E + 04 

3.00E + 02 

31 .90% 

0.012 

10.80 


1.00E + 03 

34.10% 

0.029 

18.50 


3.00E + 03 

43.20% 

0.070 

30.90 


1.00E + 04 

41.10% 

0.1 15 

47.90 

3.00E + 04 

3.00E + 02 

28.50% 

0.005 

14.90 


1.00E+03 

26.80% 

0.008 

26.40 


3.00E+03 

25.00% 

0.016 

45.30 


1.00E+04 

24.20% 

0.041 

79.50 

1.00E + 05 

3.00E + 02 

31 .70% 

0.005 

14.60 


1.00E + 03 

23.30% 

0.006 

27.80 


3.00E + 03 

20.20% 

0.008 

46.00 


1.00E+04 

19.30% 

0.020 

75.00 

3.00E + 05 

3.00E+03 

22.50% 

0.016 

10.20 


15 




Doing the Calculations 


■ A work in progress 

■ Two methods 

- LISA Calculator 

- Handy dandy web-based interface, with user selectable parameters 

- Doesn’t include spin effects 

- Calculates one source scenario 

- Problems with large (>100) mass ratios 

- Ryan Lang & Scott Hughes, does include spin effects 

- Laborious 

- Includes spin effects 

- Monte Carlo approach to spin, sky orientation, merger time 

- Problems with equal mass 

• Not included 

- Merger and ring down 

- Comprehensive exploration 



LISA Calculator 



Within the time of observation chosen, the binary system would be too relativistic for the 2PN assumptions 
used by the calculator. The observation time was changed to 0.999982 years. 


The source is chirping. The uncertainties in the set of nine parameters are: 


Parameter 

Parameter Value 

Parameter Uncertainty 

Latitude (9) 

-20° 

26.1531° 

Azimuthal Location (<p) 

305° 

45.741' 

Coalescence Time (tj 

1 year 

9.88254c 006 year 

Luminosity Distance (D L ) 

1 .06255c+008 kpc 

635608c+007 kpc 

Inclination (t) 

90° 

17.182' 

Polarization Angle (y) 

0 q 

13.1357' 

Chirp Mass (M c ) 

4607.7 M sun 

0320249 Msun 

Reduced Mass (p) 

2307 .69 

5.97818 M sim 

Orbital Phase (<p 0 ) 

0* 

171.618“ 

Sky Angle (sqrt(AO)) 

— 

593' 


Signal to noise ratio of the source: 28.6833 


The following information is included for reference: 

Observation Time: 0.999982 years 

Initial Frequency: 0.268265 mHz 

Final Frequency: 16.0827 mHz 

Frequency change (in units of f m « l/(year)); 499075 

Time to coalescence: I years 

Luminosity Distance: 1 06255e+008 kpc 

Rcdshift: 10 



Lang & Hughes 


■ Lang & Hughes calculation 

- Full 2 PN waveform 
simulation 

- Sky and polarization 

averaged 

- 1 and 2 interferometers 

- Monte Carlo spins 

- Median performance 
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Performance, using 


the LISA Calculator 


Ml 

M2 

D L 

1 IFO 

D_L D_L 

Uncertainty Uncertainty 

1 IFO 
SNR 

Initial 

Frequency 

(mHz) 

Final 

Frequency 
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0.268 

16.083 
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Distance Uncertainty 
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Lessons learned (so far) 


• Luminosity distance is the most fragile parameter. 

• Distance uncertainty doesn’t track SNR! 

• Mass is a phase measurement 

• Distance is an amplitude measurement. 

■ Inclusion of spins gives factors of 1-6 improvement. 

■ Inclusion of merger and ring-down phases may give an 
additional factor of 3. 

■ We haven’t included the hierarchical merger trees yet. 
Including spins. 



Summary 


■ Converted the science requirements from SNR based 
detection to parameter estimates from observations 

■ The adequacy of the instrument performance can only be 
verified by the forward calculation showing that the science 
requirements can be met. The instrument performance 
cannot be directly calculated from the science requirements. 

■ The performance of the instrument and the discovery space 
can be represented by “contours” of parameter estimation in 
an m 1 /m 2 coordinate system. 
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